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Basis  Construction  for  the  Design  of  Boundary  Feedback  Controls 

from  Reduced  Order  Models 

R.  Chris  Camphouse  1 
Control  Design  and  Analysis  Branch 
Air  Vehicles  Directorate 
Wright- Patters  on  Air  Force  Base 


Abstract 

In  this  paper,  we  develop  a  reduced  basis  construction  method  that  allows  for  separate  consideration  of  baseline  and  actuated 
dynamics  in  the  reduced  modeling  process.  A  prototype  initial  boundary  value  problem,  governed  by  the  two-dimensional 
Burgers  equation,  is  formulated  to  demonstrate  the  utility  of  the  method  in  a  boundary  control  setting.  A  weak  formulation 
approach,  in  combination  with  proper  orthogonal  decomposition  and  Galerkin  projection,  is  used  to  develop  a  reduced  model  of 
the  distributed  parameter  system.  Comparisons  are  done  between  reduced  and  full  order  solutions  under  open-loop  boundary 
actuation  to  illustrate  advantages  gained  by  separate  consideration  of  actuated  dynamics  in  the  reduced  modeling  process.  A 
tracking  control  problem  is  specified  for  the  full  order  system  using  a  linear  quadratic  regulator  formulation.  Comparisons  of 
feedback  control  effectiveness  are  done  to  demonstrate  benefits  in  control  effectiveness  obtained  from  separate  consideration  of 
actuated  dynamics  during  model  reduction. 
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Mb,  Ma 

tv 
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A,  v 
A,B 
G,F 

Q,R 
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covariance  matrix 

solution  snapshot 

total  number  of  snapshots 

actuated  snapshot 

total  number  of  actuated  snapshots 

basis  mode  and  temporal  coefficient,  respectively 

total  number  of  basis  modes 

total  number  of  baseline  and  actuator  modes,  respectively 

baseline  and  actuator  mode,  respectively 

baseline  and  actuator  mode  temporal  coefficient,  respectively 

eigenvalue  and  eigenvector,  respectively 

state  and  control  matrices,  respectively 

nonlinear  and  forcing  matrices,  respectively 

state  and  control  weight  matrices,  respectively 

feedback  gain  matrix 

control  robustness  parameter 

spatial  domain 

time 

spatial  coordinate  vector 
spatial  step-size 


Introduction 


Reduced  order  modeling  has  received  significant  research  attention  in  recent  years.  For  many  problems 
of  practical  interest,  the  order  of  the  system  describing  the  application  must  be  reduced.  An  illustrative 
example  where  this  is  required  is  the  development  of  feedback  control  laws  for  fluid  flow  configurations.  It 
is  not  uncommon  for  discretized  flow  models  to  describe  millions  of  state  variables.1  Unfortunately,  the 
development  of  systematic  feedback  control  laws  from  systems  of  such  large  dimension  is  a  computationally 
intractable  problem.  For  example,  if  one  uses  a  linear  quadratic  regulator  (LQR)  control  formulation,  roughly 
1012  Riccati  unknowns  need  to  be  calculated  for  a  discretized  flow  model  describing  106  states.  The  Riccati 
unknowns  are  solutions  to  a  nonlinear  matrix  equation.2  Existing  computing  power  and  computational 
algorithms  are  not  capable  of  solving  an  LQR  problem  of  such  large  dimension.  For  dynamical  models  that 

1  Mathematician,  Member  AIAA 
2210  Eighth  St.,  Bldg  146,  Room  305 
WPAFB  OH,  45433 
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are  very  large  scale,  such  as  those  describing  fluid  flow  configurations,  it  is  apparent  that  the  order  of  the 
system  must  be  reduced  prior  to  control  law  design.3 

Several  research  efforts  have  been  concerned  with,  or  relied  upon,  the  development  of  order  reduction 
strategies  that  provide  reduced  order  models  in  a  form  amenable  to  state-space  feedback  control  law  de¬ 
sign.4-14  For  the  case  of  boundary  control,  the  development  of  such  models  has  been  an  open  problem. 
For  many  applications  of  practical  interest,  such  as  feedback  control  of  the  air  flow  over  an  airplane  wing, 
boundary  actuation  is  a  requirement.  These  applications  require  actuation  to  be  located  on  the  surface  if 
they  are  to  be  implemented  in  hardware  in  the  physical  system.  Most  reduced  modeling  efforts  have  either 
been  concerned  with  the  case  of  control  action  via  a  body  force  or  have  approximated  boundary  actuation 
by  a  body  force  on  the  domain  interior.  The  reduced  modeling  task  is  usually  simpler  in  the  case  of  body 
forcing  due  to  the  presence  of  homogeneous  boundary  conditions. 

Few  techniques  have  been  obtained  thus  far  providing  reduced  order  models  with  explicit  boundary 
input  that  can  be  used  for  systematic  design  of  feedback  control  laws.  The  case  of  boundary  actuation  is 
simply  more  difficult  than  the  case  of  body  forcing.  It  is  difficult  to  extract  the  boundary  control  input  in 
the  reduced  model.  When  proper  orthogonal  decomposition  (POD)  is  used  in  conjunction  with  Galerkin 
projection,  boundary  conditions  are  absorbed  in  the  process.  An  additional  step  is  required  to  make  the 
action  of  the  control  input  explicit  in  the  model.  Approximating  boundary  actuation  as  a  body  force  near  the 
boundary  is  convenient.  It  circumvents  the  difficulty  associated  with  making  boundary  control  inputs  explicit 
in  the  reduced  model.  However,  feedback  control  by  body  forcing  is  a  very  different  problem  mathematically 
than  the  problem  of  boundary  control.  In  the  context  of  systems  governed  by  partial  differential  equations, 
boundary  control  leads  to  an  unbounded  control  input  operator.15  In  other  words,  when  viewed  as  an 
operator  equation,  boundary  control  leads  to  an  unbounded  B  operator  in  the  state  equation 

x  =  Ax  +  Bu ,  (1) 

x(0)  =  x0.  (2) 

For  the  case  of  body  forcing,  the  B  operator  is  typically  bounded.  As  a  result,  while  convenient,  the  reliability 
associated  with  approximating  a  boundary  control  as  a  distributed  body  force  is  questionable. 

Recently,  inroads  have  been  made  allowing  for  the  extraction  of  boundary  conditions  in  reduced  order 
POD  models.  For  the  most  part,  these  methods  utilize  POD  in  combination  with  a  weak  formulation  of  the 
Galerkin  projection.16  19  The  advantage  of  a  weak  formulation  is  that  boundary  conditions  appear  naturally 
when  the  reduced  model  is  written  weakly.  In  this  paper,  we  utilize  a  method16  incorporating  difference 
approximations  in  the  weak  formulation  to  construct  reduced  models  with  boundary  control  input  appearing 
explicitly. 

Special  care  must  be  taken  when  POD  is  used  for  the  development  of  feedback  control  laws.  For  many 
canonical  control  formulations,  the  dynamics  induced  by  the  feedback  control  are  not  known  a  priori.  Thus, 
it  is  difficult  to  construct  a  set  of  POD  basis  modes  capable  of  spanning  the  baseline  solution  as  well  as  the 
dynamics  of  the  system  with  boundary  actuation.  It  is  entirely  possible  that  the  dynamics  induced  by  a 
feedback  control  are  not  adequately  represented  by  the  reduced  model  if  the  POD  basis  is  not  constructed 
carefully.  Construction  of  a  POD  basis  capable  of  spanning  the  baseline  solution,  as  well  as  dynamics 
introduced  by  boundary  actuation,  is  the  subject  of  this  paper. 

Model  Problem 


A  distributed  parameter  system  is  formulated  which  models  convective  flow  over  an  obstacle.  The  re¬ 
sulting  model  problem  is  used  to  demonstrate  the  importance  of  careful  reduced  model  design.  Let  Hi  CM2 
be  the  rectangle  given  by  (a,  b]  x  (c,  d).  Let  D2  C  Di  be  the  rectangle  given  by  [ai,a2]  x  [&i ,  62]  where 
a  <  a\  <  a2  <  b  and  c  <  b\  <  62  <  d.  The  problem  domain,  D,  is  given  by  D  =  Di  \  D2.  In  this  configura¬ 
tion,  D2  is  the  obstacle.  Dirichlet  boundary  controls  are  located  on  the  obstacle  bottom  and  top,  denoted 
by  T#  and  T^,  respectively. 

The  dynamics  of  the  system  are  described  by  the  two-dimensional  Burgers  equation 

d  1 

—w(t,x,y)  +  V  •  F(w)  =  —  Aw(t,x,y)  (3) 
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Figure  1:  Problem  Geometry. 


for  t  >  0  and  (x,y)  G  fh  In  (3),  F(w)  has  the  form 

F(w)  = 

where  C±,  C2  are  nonnegative  constants.  This  equation  has  a  convective  nonlinearity  like  that  found  in 
the  Navier-Stokes  momentum  equation  modeling  fluid  flow.20  The  quantity  Re,  a  nonnegative  constant,  is 
analogous  to  the  Reynolds  number  in  the  Navier-Stokes  momentum  equation. 

To  complete  the  model  of  the  system,  boundary  conditions  must  be  specified  as  well  as  an  initial  condition. 
For  simplicity,  boundary  controls  are  assumed  to  be  separable.  With  this  assumption,  we  specify  conditions 
on  the  obstacle  bottom  and  top  of  the  form 

w(t,VB)  =  uB(t)^B  (x),  (5) 

w(t,rT)  =  uT(t)^T{x).  (6) 

In  (5)- (6),  uB(t)  and  ur(t)  are  the  controls  on  the  bottom  and  top  of  the  obstacle,  respectively.  The  profile 

functions  ^B(x)  and  Tt(x)  describe  the  spatial  influence  of  the  controls  on  the  boundary.  A  parabolic  inflow 
condition  is  specified  of  the  form 

w(t,Tin)  =  f(y).  (7) 

At  the  outflow,  a  Neumann  condition  is  specified  according  to 

Tout)  =  0.  (8) 

For  notational  convenience,  denote  the  remaining  boundary  as  Tjy.  We  require  that  values  be  fixed  at  zero 
along  Tv  as  time  evolves.  The  resulting  boundary  condition  is  of  the  form 

w(t,Tu)  =  0 

The  initial  condition  of  the  system  is  given  by 

w(0,x,y)  =  w0(x,y)  e  L2(n).  (10) 


C\ 


C2 


wz(t,x,y) 


(4) 


Reduced  Basis  Construction 


In  this  work,  we  use  a  POD  algorithm21  based  on  the  snapshot  method22  to  construct  the  low  order 
basis  needed  for  the  development  of  the  reduced  order  model.  A  data  ensemble  of  snapshots  is 

generated  for  the  system  via  numerical  simulation  or  experiment,  where  N  is  the  total  number  of  snapshots. 
Each  snapshot  consists  of  instantaneous  system  data.  With  the  snapshot  ensemble  in  hand,  the  N  x  N 
correlation  matrix  L  defined  by 

L^^iSuSj)  (11) 

is  constructed.  In  this  work,  we  utilize  the  standard  L2(Q)  inner  product 

[  SiS; dx,  (12) 

Jn 
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where  S*  denotes  the  complex  conjugate  of  Sj,  in  the  construction  of  L. 

The  eigenvalues  {A^}^  of  L  are  calculated  and  sorted  in  descending  order.  The  ratio 


100 


(z£ i 
VET 


A,: 


(13) 


is  used  to  determine  the  number  of  POD  basis  functions  to  construct.  The  quantity  in  (13)  provides  a 
measure  of  the  ensemble  energy  that  is  captured  by  a  POD  basis  consisting  of  M  modes.  By  requiring  a 
percentage  of  the  energy  contained  in  the  snapshot  ensemble  be  contained  in  the  basis,  the  smallest  value  of 
M  is  calculated  such  that  the  quantity  in  (13)  is  greater  than  or  equal  to  that  percentage. 

The  eigenvectors  corresponding  to  the  M  eigenvalues  of  largest  magnitude  are  calculated.  Each 

eigenvector  is  normalized  so  that 

IM2  =  y.-  (14) 

The  orthonormal  POD  basis  set  {0i(x)}^1  is  constructed  according  to 


N 

<MX)  = 

3= 1 


(15) 


where  Vij  is  the  jth  component  of  Vi.  The  POD  basis  is  optimal  in  an  energy  sense.21  It  captures  the  mean 
square  energy  of  the  snapshot  ensemble  better  than  any  other  basis. 

With  the  basis  in  hand,  the  system  solution  re(t,x)  is  approximated  as  a  linear  combination  of  POD 
modes,  i.e., 

M 

w(t,x.)  «  y^ai(*)0i(x).  (16) 

i=  1 

Galerkin  projection  onto  the  basis  results  in  a  system  of  ordinary  differential  equations  for  the  temporal 
coefficients 


Snapshot  Decomposition 


Using  an  energy  argument  based  on  (13)  is  a  convenient  way  to  determine  the  number  of  POD  modes 
needed  for  the  reduced  model.  However,  blindly  applying  (13)  to  determine  the  necessary  number  of  modes 
is  problematic  in  a  feedback  control  setting.  First,  the  snapshot  ensemble  typically  consists  of  baseline  (zero 
control  input)  and  actuated  data.  Actuated  data  is  obtained  by  stimulating  the  system  with  open-loop 
control  inputs.  Using  an  energy  argument  based  on  open- loop  information  to  make  conclusions  about  the 
energy  of  the  system  under  feedback  control  is  problematic.  The  energy  content  of  open- loop  data  is  likely 
very  different  than  the  energy  content  of  the  system  under  feedback  control.  Second,  in  many  applications  it 
is  possible  that  energy  content  in  the  snapshot  ensemble  is  dominated  by  that  in  the  baseline  data.  This  is 
particularly  troublesome  in  the  context  of  boundary  control.  In  boundary  control  applications,  it  is  usually 
desired  that  control  input  energy  be  as  small  as  possible  while  still  satisfying  the  control  objective.  For 
example,  in  flow  control  applications  where  the  control  is  located  on  a  surface,  it  is  desired  that  small  control 
inputs  yield  large  changes  in  the  flow  field  behavior.23  In  essence,  effects  of  small  control  inputs  are  amplified 
by  the  natural  instabilities  and  comparatively  high  energy  content  of  the  baseline  flow  field.  Simply  applying 
(13)  to  an  ensemble  consisting  of  baseline  and  actuated  data  presents  the  risk  of  important  structures  due 
to  control  input  being  discarded  if  they  are  of  much  lower  energy  content  than  structures  associated  with 
the  baseline  solution.  For  these  reasons,  we  extend  the  POD  algorithm  described  above  so  that  baseline  and 
control  input  energy  are  considered  separately.  The  resulting  basis  set  will  consist  of  modes  significant  to 
the  baseline  solution  as  well  as  those  significant  to  the  system  under  control  actuation. 

The  basic  idea  is  to  decompose  each  snapshot  in  the  ensemble  into  a  component  in  the  span  of  the 
baseline  POD  basis  and  an  orthogonal  component.  This  is  done  by  employing  useful  properties  of  orthogonal 
projections  on  Hilbert  spaces.24  By  considering  the  case  of  baseline  and  actuated  data  separately,  the 
orthogonal  component  is  constructed  such  that  it  contains  new  information  due  to  the  control  input.  This 
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new  information  can  be  used  in  conjunction  with  (13)  to  determine  additional  modes  that  are  significant 
from  a  control  standpoint. 

An  ensemble  of  solution  snapshots  is  generated  for  the  case  of  zero  control  input.  From  this  baseline 
snapshot  ensemble,  a  set  of  POD  modes  is  constructed  using  the  snapshot  algorithm  discussed  previously. 
For  notational  convenience,  denote  the  baseline  basis  by  where  Mb  is  the  total  number  of  baseline 

modes.  By  employing  the  energy  ratio  in  (13),  Mb  can  be  chosen  such  that  the  baseline  basis  contains  an 
arbitrary  amount  of  the  energy  contained  in  the  baseline  snapshot  ensemble. 

With  a  set  of  baseline  modes  in  hand,  an  ensemble  of  solution  snapshots  is  generated  for  the  case  of 
nonzero  control  input.  Denote  the  actuated  ensemble  by  {Ti]Mx  where  Na  is  the  total  number  of  actuated 
snapshots. 

For  each  snapshot  in  the  actuated  ensemble,  we  determine  the  component  that  is  in  the  span  of  the 
baseline  basis.  In  particular,  define  b^  according  to 

bij  =  {Ti,£j) ,  1  <  i  <  Na ,  1  <  j  <  MB .  (17) 

Then,  b^  is  the  projection  of  the  ith  actuated  snapshot  Ti  onto  the  j th  baseline  POD  mode  In  other 
words,  the  product  bij^j  is  the  component  of  T{  that  is  in  the  direction  of  £j.  The  linear  combination 
Yl!f=i  bij€j  is  the  component  of  Ti  in  the  span  of  the  baseline  basis. 

Define  T)  according  to 

Mb 

Ti  =  Ti  —  bij£j.  (18) 

3  = 1 

Then,  Ti  is  the  component  of  Ti  not  contained  in  the  span  of  the  baseline  basis.  As  Ti  is  a  solution  snapshot 
for  the  case  of  nonzero  control  input,  %  consists  of  new  information  due  to  the  control  input  that  is  not 
already  contained  in  the  span  of  the  baseline  modes. 

A  second  set  of  POD  modes  is  constructed  from  the  data  ensemble  Denote  this  set  of  “actuator 

modes”  by  where  Ma  is  the  total  number  of  modes.  The  energy  ratio  in  (13)  is  used  to  determine 

Ma  such  that  the  basis  of  actuator  modes  contains  an  arbitrary  amount  of  the  additional  energy  resulting 
from  the  control  input. 

By  construction,  the  set  of  baseline  modes  consists  of  orthonormal  basis  functions.  Similarly,  the  actuator 
modes  are  orthonormal.  It  is  easy  to  show  that  the  set  of  baseline  modes  is  orthogonal  to  the  set  of  actuator 
modes.  Consider  the  inner  product  of  £&  and  Ti  for  arbitrary  i,k.  We  see  that 

/mb  \ 

=  (/Ti  -  y ~'/ij€j,€k) 

Mb 

J  =  1 

As  the  baseline  modes  are  orthonormal,  (£j,£fc)  =  0  unless  j  =  k.  For  j  =  fc,  (£j,£k)  =  1.  As  a  result, 


(19) 

(20) 


Mb 

{Ti,tik)  -  bij  (£j,€k)  =  (Ti,£k)  -  bik •  (21) 

3  = 1 

By  definition  (17),  b^  =  (T^,^).  Thus,  we  have 

(Ti^k^  =  0,  for  arbitrary  i,  k.  (22) 

From  (15),  each  actuator  mode  rji  is  a  linear  combination  of  the  snapshots  {T^}^.  As  a  result,  (22)  in 
combination  with  the  linearity  of  the  inner  product  yields 

{ViiCk}  =  for  arbitrary  i,k.  (23) 
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Therefore,  we  have 


(24) 


{Vi}  -L  {0 }• 

This  result  allows  us  to  combine  the  baseline  and  actuator  modes  into  an  overall  basis  set 


(25) 


where  all  modes  in  the  basis  are  orthonormal.  This  is  advantageous  for  the  development  of  reduced  order 
models.  The  system  solution  w(t,x)  is  still  approximated  as  a  linear  combination  of  modes  as  in  (16). 
Moreover,  separate  consideration  of  baseline  and  actuated  energy  allows  us  to  write  this  linear  combination 
as 

Mb  Ma 

w(t,x)  «  (*)&(*)  +  5Z&W»?i(x),  (26) 

3= 1  'i'=l 


where  and  are  temporal  coefficients  for  the  baseline  and  actuator  basis,  respectively.  This 

allows  us  to  consider  the  system  solution  as  a  baseline  component  and  an  additional  component  induced  by 
the  control  input. 

Reduced  Order  Model 

We  now  develop  a  reduced  order  model  for  the  system  described  by  (3),  (5)- (10).  Our  aim  is  to  develop 
a  low  order  state-space  model  where  the  Dirichlet  control  inputs  specified  by  (5)- (6)  appear  in  the  model 
explicitly. 

Taking  the  inner  product  of  both  sides  of  (3)  with  the  ith  POD  mode  <fii(x,y )  and  utilizing  Green’s 
identities  results  in  the  weak  formulation 


/ 

Jn 


d  1 

—w(t,x,y)<f>i(x,y)dx  =  — 


IdQ 


dn 


(Vw(t,x,y)  ■  n)(f>i(x,y)dA(x)  -  /  Vw(t,x,y)  ■’V4>i(x,y) d> 

Jn 

(F(w)  ■  n)(f>i(x,y)dA(x)  -  [  F-V<f>i(x,y) d> 

Jn 


(27) 


where  n  denotes  the  unit  outward  normal. 

As  seen  in  (15),  each  POD  mode  is  a  linear  combination  of  solution  snapshots.  From  (9),  snapshot  values 
along  Yu  are  specified  to  be  zero.  As  a  result,  POD  modes  are  zero  along  Yu-  Thus,  the  first  boundary 
integral  in  (27)  is  decomposed  as 


dn 


<yw(t,  x,  y)-n)4>i(x,  y)dA(x) 


nd  d  \  fd  d 

— w(t,x,bi)(t>i(x,bi)  -  —w(t,x,b2)Mx,b2)J  dx  -  J  — w(t,a,y)<j>i(a,y)dy , 


(28) 


where  condition  (8)  has  been  used  to  specify  that 

fd  d 

J  ~dxW^' b ’  y^dy  = 

In  a  similar  fashion,  the  remaining  boundary  integral  in  (27)  is  decomposed  as 


(29) 


/  (F(w)  ■  n)<f)i(x,y)dA(x)  =  -  /  (w(t,b,y)2(t>i(b,y)  -  f{y)2<t>i{a,y))  dy, 

Jdn  z  Jc 


(30) 


where  (7)  has  been  used  to  incorporate  the  inflow  condition  f(y). 

Control  inputs  and  the  Dirichlet  inflow  condition  are  not  explicit  in  (28).  They  can  be  made  explicit  by 
approximating  partial  derivatives  along  the  boundary.  For  h  >  0,  we  see  that 


9  uB{t)^B{x) -w{t,x,bi  -  h) 

—w(t,x,b i)  « - - - , 

dy  h 

d  w(t,x,b2  +  h)  -ur(t)yT(x) 

—w(t,x,b2)  « - 7 , 

dy  h 

d  „  tU  „  „  x  _  w(t,  a  +  h,y)  —  f(y) 

— w(t,  a,  y)  w - -h - . 


(31) 

(32) 

(33) 
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These  expressions  are  substituted  into  (28).  Approximating  w(t,x,y)  as  a  linear  combination  of  POD 
modes  in  (27),  (28),  and  (30)  results  in  a  reduced  order  system  model.  This  system  is  of  the  form 

d  =  Ad  T  Bu  T  G{d )  T  T1,  (34) 


where 


1  fa2 

A(i,j)  =  J  -  h)cj)i(x,bi)  +  (j)j(x,b2  +  h)</>i(x,  b2))dx 


+ 


/  ct)j(a  +  h1y)(l)i(a,y)dy  +  h  /  V^(x, y)  •  V<^(z, 2/)]dx, 
Jc  J  O 


(35) 


B  = 


hRe 


/“2  0i (a:,  &i)^B(a:)da;  /“2  0i(x,  M’MaOcte 


LC  4>m(x,  bi)^B{x)dx  /“2  4>m{x,  b2)^T{x)dx 


Mx  2 


(36) 


G(a)  =  5 


In  (j2jLi  aj Mx^y))  ^0i(a;,y)dx  -  /Cd  (Ejli  ®j<Pj(.b,y ))  <t>i{b,y)dy 


In  (Eji i  oij<t>j{x, y))  £<Pm(x, y) dx  -  /cd  (E^li  aj4>j(b , y))  <M&,  y)dy 


(37) 


Mx  1 


r 


F  = 


/c  (hkf(y)  +  y(y)2)Ma,y)dy 


\Jc(hkf(y)  +  |/(2/)2)0M(a,y)dyJ 


Mx  1 


(38) 


= 


^s(t) 

ur(t) 


2x1 


(39) 


Projecting  the  initial  condition  reo(x,  y)  onto  the  POD  basis  results  in  an  initial  condition  for  the  reduced 
order  model  of  the  form 

<a(0)  =  dQ.  (40) 


Open-Loop  Comparisons 


We  now  demonstrate  the  impact  of  basis  construction  technique  on  the  ability  of  the  reduced  model  to 
represent  dynamics  induced  by  a  boundary  control  input.  Instantaneous  snapshots  are  generated  for  (3), 
(5)- (10)  via  numerical  simulation.  A  positive  parabolic  profile  with  unit  maximum  amplitude  is  specified  for 
the  inlet  condition  in  boundary  condition  (7).  In  (4),  we  set  C i  =  1  and  C2  =  0  in  order  to  obtain  solutions 
that  convect  from  left  to  right  for  the  positive  inlet.  In  addition,  we  specify  that  Re  =  300.  The  problem 
domain  Q  is  discretized,  resulting  in  a  uniform  grid  with  spatial  step-size  h.  We  utilize  a  finite  difference 
scheme25  to  numerically  solve  the  model  problem  with  and  without  boundary  control  input.  The  resulting 
discretized  system  describes  roughly  2,000  states. 
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We  construct  reduced  order  models  from  snapshot  ensembles  obtained  for  two  scenarios.  In  the  first 
scenario,  snapshots  are  generated  for  the  baseline  solution  and  for  boundary  actuation  at  a  fixed  frequency. 
For  the  second  case,  snapshots  are  generated  for  the  baseline  solution  and  a  more  complicated  boundary 
input  where  the  frequency  varies  with  time.  For  both  scenarios,  we  compare  model  agreement  resulting  from 
combining  baseline  and  actuator  snapshots  into  an  overall  lumped  snapshot  ensemble  to  that  obtained  by 
decomposing  snapshots  into  their  baseline  and  actuated  components  and  constructing  the  POD  basis  as  in 
(25).  In  the  results  that  follow,  basis  construction  from  a  lumped  snapshot  ensemble  is  referred  to  as  the 
“lumped”  method.  Constructing  the  basis  by  considering  baseline  and  actuated  energy  separately  is  referred 
to  as  the  “split”  method. 

Scenario  1 


Snapshots  are  generated  for  the  baseline  solution  and  for  the  solution  arising  under  periodic  boundary 
actuation.  Inputs  specified  are  of  the  form 

UB(t)  =  sin(irt)  ur(t)  =  0,  (41) 

u  b  (£)  =  0  uT(t)  =  sin{nt).  (42) 

In  (41),  periodic  actuation  is  done  on  the  bottom  of  the  obstacle  while  values  along  the  obstacle  top  are  held 
at  zero.  In  (42),  values  along  the  obstacle  bottom  are  held  at  zero  with  periodic  actuation  occurring  at  the 
top.  For  each  control  input  listed  in  (41)- (42),  snapshots  are  taken  in  increments  of  At  =  0.1  starting  from 
£  it  0  and  ending  at  £  =  15.  The  steady  baseline  solution  is  used  for  the  initial  condition. 

With  an  ensemble  of  snapshots  in  hand,  ratio  (13)  is  used  to  determine  the  number  of  basis  modes  to 
construct.  Requiring  that  99.9%  of  the  ensemble  energy  be  contained  in  the  POD  basis  results  in  a  lumped 
basis  consisting  of  7  modes.  Separate  consideration  of  baseline  and  actuated  energy  results  in  a  split  basis 
consisting  of  1  baseline  mode  and  16  actuator  modes. 

We  now  employ  linear  combination  (16)  to  compare  boundary  condition  agreement  between  the  full  order 
system  and  reduced  models  obtained  via  the  lumped  and  split  POD  bases.  By  specifying  characteristic 
functions  for  the  control  profiles  ^(x)  and  Tt(^)  in  (5)-(6),  we  see  that 

M 

^  (£)</>*  (T# )  ~  ?/;(£,  T^)  =  ub (£) 5 

i= 1 
M 

^  ai(t)(j)i(YT )  ~  w{t,  Tt)  =  uT{t ), 

i= 1 

We  construct  the  linear  combinations  on  the  left  in  (43)- (44)  and  compare  them  to  the  exact  boundary 
conditions  us{t)  and  ur(t)  specified  in  the  full  order  system.  We  first  perform  comparisons  for  boundary 
inputs  explicitly  used  during  ensemble  generation.  Results  obtained  for  the  baseline  solution  and  for  the 
solution  with  periodic  boundary  actuation  of  the  form  sin{i:t )  are  shown  in  Figure  2.  In  that  figure,  dashed 
curves  denote  the  linear  combination  of  POD  modes  restricted  to  the  boundary.  Solid  curves  denote  the 
exact  full  order  boundary  input.  Results  obtained  for  the  lumped  basis  method  are  plotted  on  the  left. 
Split  POD  basis  results  are  plotted  on  the  right.  As  seen  in  Figure  2,  both  methods  result  in  very  good 
agreement  between  the  exact  boundary  conditions  and  the  linear  combination  of  POD  modes  restricted  to 
the  boundary.  Dashed  and  solid  curves  are  virtually  identical. 

In  a  feedback  control  setting,  dynamics  induced  by  the  control  are  typically  not  known  a  priori.  It  is  not 
possible  to  specify  a  boundary  feedback  control  law  as  an  open-loop  input  during  ensemble  creation.  The 
specifics  of  the  control  law  and  the  dynamics  induced  by  it  are  not  known  at  that  stage.  Typically,  in  the 
closed-loop  system,  the  boundary  input  resulting  from  the  control  law  will  be  different  than  the  inputs  used  to 
generate  the  snapshot  ensemble.  As  a  result,  it  is  useful  to  compare  reduced  and  full  order  model  agreement 
for  inputs  not  specified  during  ensemble  creation.  This  provides  insight  into  the  suitability  of  the  reduced 
model  for  closed-loop  control  law  design.  For  these  reasons,  we  now  compare  boundary  condition  agreement 
between  the  reduced  and  full  order  systems  for  open-loop  inputs  that  were  not  used  during  ensemble  creation. 


(43) 

(44) 
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Figure  2:  Boundary  Condition  Accuracy. 


Boundary  inputs  specified  are  of  the  form 


uB{t)  =  min  1  j  , 

ur(t)  =  sin  ^71 't^j  . 

Results  obtained  for  the  lumped  and  split  basis  methods  are  shown  in  Figure  3. 

Lumped  Split 


(45) 

(46) 


Figure  3:  Boundary  Condition  Accuracy. 

As  seen  in  Figure  3,  very  good  agreement  is  seen  between  the  reduced  and  full  order  systems  for  the 
split  method,  even  though  the  inputs  considered  were  not  specifically  included  in  the  snapshot  ensemble. 
Condition  (46)  is  reconstructed  well  using  the  lumped  method.  However,  the  reconstruction  of  the  piecewise 
linear  input  in  (45)  is  much  less  accurate  when  the  lumped  basis  is  used. 

To  further  compare  the  lumped  and  split  basis  methods  and  their  utility  for  control  law  design,  we  project 
the  full  order  solution  at  each  time  step  onto  the  lumped  and  split  POD  bases.  The  resulting  temporal 
coefficients  are  compared  to  those  predicted  by  the  reduced  order  models.  Boundary  inputs  specified  are 
as  in  (45)- (46).  Results  obtained  for  the  first  5  temporal  coefficients  of  the  lumped  method  are  shown  in 
Figure  4.  The  first  5  temporal  coefficients  for  the  split  method  are  shown  in  Figure  5.  In  Figures  4-5, 
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Figure  4:  Temporal  Coefficient  Accuracy  for  the  Lumped  Method. 

solid  curves  denote  values  of  temporal  coefficients  obtained  via  the  projection.  Dashed  curves  denote  the 
solution  of  the  reduced  order  model.  Overall,  the  accuracy  of  the  split  method  is  better,  particularly  for 
temporal  coefficients  with  higher  frequency  content.  Separate  consideration  of  actuated  energy  during  basis 
construction  results  in  better  representation  of  dynamics  induced  by  boundary  inputs  not  specified  during 
ensemble  creation.  This  is  advantageous  in  a  boundary  feedback  control  setting  where  dynamics  induced  by 
the  control  are  typically  not  known  beforehand. 


Split 


Figure  5:  Temporal  Coefficient  Accuracy  for  the  Split  Method. 


Scenario  2 

It  is  likely  that  a  snapshot  ensemble  for  inputs  at  a  single  frequency  results  in  a  POD  basis  that  does 
not  adequately  span  the  dynamics  induced  by  a  feedback  control.  Feedback  controls  designed  from  such  a 
basis  are  bound  to  be  ineffective  when  implemented  in  the  full  order  system.  In  Scenario  2,  we  compare  the 
lumped  and  split  basis  methods  using  a  snapshot  ensemble  generated  from  boundary  inputs  that  excite  a 
range  of  system  dynamics.  Inputs  specified  are  of  the  form 

usit)  =  sin  ^7 r  (e*)0’3^  ur{t)  =  0, 

u b (t)  =  0  uT(t)  =  sin  (it  (et)°‘3 
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(47) 

(48) 


The  sinusoidal  input  function  is  shown  in  Figure  6.  As  seen  in  that  figure,  an  input  of  this  form  generates 
the  system  response  over  a  range  of  frequencies.  The  resulting  POD  basis  is  much  more  likely  to  adequately 
span  the  unknown  dynamics  generated  by  a  feedback  control. 

Instantaneous  snapshots  are  generated  for  the  baseline  solution  as  well  as  for  solutions  arising  from  the 
inputs  in  (47)- (48).  For  each  control  input  listed  in  (47)- (48),  snapshots  are  taken  in  increments  of  At  =  0.1 
starting  from  t  —  0  and  ending  at  t  =  15.  The  steady  baseline  solution  is  used  for  the  initial  condition. 

Lumped  Split 


Figure  7:  Boundary  Condition  Accuracy. 

Requiring  that  99.9%  of  the  ensemble  energy  be  contained  in  the  POD  basis  results  in  a  lumped  basis 
consisting  of  7  modes.  A  split  basis  comprised  of  1  baseline  mode  and  25  actuator  modes  are  needed  when 
baseline  and  actuated  energy  is  considered  separately. 

As  in  Scenario  1,  we  compare  boundary  condition  accuracy  of  the  lumped  and  split  methods  for  boundary 
inputs  not  specified  during  ensemble  generation.  For  the  sake  of  comparison,  we  use  the  boundary  conditions 
given  by  (45)- (46).  Results  obtained  for  the  lumped  and  split  basis  methods  are  shown  in  Figure  7.  As  seen 
in  that  figure,  very  good  boundary  condition  agreement  is  seen  between  the  reduced  and  full  order  systems 
for  both  basis  methods.  In  particular,  by  comparing  Figures  3  and  7,  we  see  that  the  piecewise  linear 
boundary  condition  in  (45)  is  represented  much  better  by  the  lumped  method  when  inputs  (47)-(48)  are 
used  to  generate  the  snapshot  ensemble. 

As  before,  we  now  project  the  full  order  solution  at  each  time  step  onto  the  lumped  and  split  POD  bases. 
The  resulting  temporal  coefficients  are  compared  to  those  predicted  by  the  reduced  order  models.  The 


11 


results  for  the  lumped  method  are  shown  in  Figure  8.  Split  method  results  are  shown  in  Figure  9.  As  seen  in 
those  figures,  the  accuracy  of  the  split  method  is  better.  For  temporal  coefficients  with  significant  frequency 
content,  values  predicted  by  the  reduced  model  are  virtually  identical  to  those  obtained  by  projecting  the  full 
solution  onto  the  split  basis,  even  though  the  boundary  inputs  specified  are  different  than  those  used  during 
ensemble  creation.  The  split  method  is  better  suited  with  regard  to  feedback  control  law  design  as  it  is  more 
capable  of  accurately  representing  dynamics  that  are  not  explicitly  included  in  the  snapshot  ensemble. 


Lumped 


Figure  8:  Temporal  Coefficient  Accuracy  for  the  Lumped  Method. 


Split 


Figure  9:  Temporal  Coefficient  Accuracy  for  the  Split  Method. 


Control  Formulation 

We  now  specify  the  control  formulation  used  to  compare  the  effectiveness  of  the  lumped  and  split  basis 

methods  in  a  feedback  control  setting.  The  reduced  system  given  by  (34),  (40)  is  linearized  about  the  origin 

yielding  a  state-space  equation  of  the  form 

a(t)  =  Aa  +  Bu ,  (49) 

a(0)  =  ctQ.  (50) 

We  consider  the  tracking  problem  for  (49)-(50).  A  fixed  reference  signal  wref{x)  is  specified  for  the 
full  order  system.  Projecting  wref{x)  onto  the  POD  basis  yields  tracking  coefficients  for  the  reduced  order 
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model,  denoted  by  aref.  The  dynamics  of  the  linearized  model  under  tracking  control  are  given  by 


a 

Oiref 


A  0 
0  0 


a 

Oiref 

AX  T  Bu , 


where  we  have  defined  the  augmented  state  X  as 


X{t)  = 


a{t) 

O ref 


Oo 

Oref 


(51) 

(52) 

(53) 


with  Xo  = 

To  formulate  the  control  problem,  we  consider  the  y-shifted  linear  quadratic  regulator  (LQR)  cost  function 

/‘OO 

J(ao,u)  =  /  {(a  —  aref)TQ(o  —  oref)  +  uT Ru]  e2ltdt.  (54) 

Jo 

In  (54),  Q  is  a  diagonal,  symmetric,  positive  semi-definite  matrix  of  state  weights.  R  is  a  diagonal,  symmetric, 
positive  definite  matrix  of  control  weights.  The  quantity  7,  a  nonnegative  constant,  is  an  additional  parameter 
that  provides  added  robustness  in  the  control.26,27  The  optimal  control  problem  we  consider  is  to  minimize 
(54)  over  all  controls  u  G  L2(0,  00)  subject  to  the  constraints  (51)-(53). 

For  a  controllable  system,  the  LQR  problem  has  a  unique  solution  of  the  form 


^ opt  K  X 

=  -\Ki  K2]X 

=  -  [r-1bth11  r-1bth12 


where  Iln  is  the  unique  symmetric,  non- negative  solution  of  the  algebraic  Riccati  equation 

(A  +  7/)Tnn  +  n U{A  +  7 1)  -  n11BBrtBTii11  +  q  =  o. 

The  matrix  II12  in  (57)  satisfies  the  equation 

[(A  +  7/)t  -  n  11BR~1Bt]  n12  =  Q. 


(55) 

(56) 

(57) 


(58) 


(59) 


The  feedback  control  obtained  from  the  linearized  model  is  placed  into  the  nonlinear  state-space  equation. 
The  resulting  closed-loop  nonlinear  system  is  of  the  form 


X  =  (A  —  BK)X  +  [G(o t)  0]T  +  [F  0]T, 

X(0)  =  x0. 


(60) 

(61) 


Closed-Loop  Results 

We  use  the  LQR  formulation  in  (54)  to  compare  closed-loop  results  obtained  from  the  lumped  and  split 
basis  methods.  A  snapshot  ensemble  is  constructed  containing  baseline  solution  data  as  well  as  data  resulting 
from  nonzero  boundary  actuation.  As  it  is  desired  that  the  POD  basis  spans  unknown  dynamics  introduced 
by  the  LQR  feedback  control,  boundary  inputs  specified  during  ensemble  creation  are  of  the  form 

ub(J)  =  Csin  ^7r  (et)°‘3^  ur(t)  =  0,  (62) 

us(t)  =  0  ur{t)  =  Csin  (^7r  (63) 


for  C=  1,2,3. 

For  each  control  input  listed  in  (62)- (63),  snapshots  are  taken  in  increments  of  At  =  0.1  starting  from 
t  —  0  and  ending  at  t  =  15.  The  steady  baseline  solution  is  specified  for  the  initial  condition.  The  resulting 
snapshot  ensemble  consists  of  roughly  900  snapshots.  Requiring  that  99%  of  the  ensemble  energy  be  contained 
in  the  POD  basis  results  in  a  lumped  basis  consisting  of  5  modes.  A  split  basis  comprised  of  1  baseline  mode 
and  20  actuator  modes  is  needed  when  baseline  and  actuated  energy  are  considered  separately.  The  first  9 
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Figure  10:  Split  Basis  Modes. 


Figure  11:  Tracking  Reference  Function. 


modes  of  the  split  basis  are  shown  in  Figure  10.  In  that  figure,  mode  1  is  the  baseline  mode.  Modes  2-9  are 
actuator  modes. 

The  tracking  LQR  problem  requires  the  specification  of  the  reference  signal  aref.  In  the  results  that 
follow,  aref  is  obtained  from  the  unactuated  steady  solution  for  the  case  Re  =  50.  This  solution  is  projected 
onto  the  lumped  and  split  bases.  The  temporal  values  obtained  are  used  as  tracking  coefficients  in  the  reduced 
order  control  problem.  The  reference  signal  obtained  by  projecting  the  steady  solution  at  Re  =  50  onto  the 
split  basis  is  shown  in  Figure  11.  The  reference  function  obtained  by  projecting  onto  the  lumped  basis  is 
similar.  To  complete  the  control  formulation,  each  state  in  the  reduced  order  model  is  prescribed  a  weight 
of  2,500.  The  two  boundary  controls  are  each  given  unit  weight.  The  value  specified  for  7  in  (54)  is  0.25. 
The  closed- loop  solution  of  the  reduced  order  model  constructed  with  the  split  POD  basis  is  shown  in  Figure 
12.  By  comparing  the  controlled  solution  of  Figure  12  to  the  reference  function  in  Figure  11,  it  is  apparent 
that  the  closed-loop  reduced  order  model  satisfies  the  control  objective  quite  well.  Separate  consideration  of 
actuated  energy  in  the  split  POD  basis  method  results  in  satisfactory  tracking  of  the  reference  signal. 

When  the  energy  ratio  in  (13)  is  applied  to  baseline  and  actuated  data  lumped  together  into  an  overall 
snapshot  ensemble,  the  results  are  much  less  favorable.  Closed-loop  solutions  of  the  reduced  order  model 
constructed  with  the  lumped  POD  basis  are  shown  in  Figure  13.  As  seen  in  that  figure,  virtually  no  tracking 
is  achieved  by  the  reduced  order  control.  Adjusting  parameters  in  the  control  formulation  has  little  effect 
on  this  result.  Increasing  the  state-weights  and  the  parameter  7  to  10,000  and  0.75,  respectively,  does  not 
significantly  improve  the  performance  of  the  control.  System  information  relevant  from  a  control  standpoint 
is  discarded  when  an  energy  argument  is  applied  during  order  reduction  to  the  lumped  snapshot  ensemble 
containing  baseline  and  actuated  data.  The  resulting  reduced  order  model  does  not  adequately  describe 
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Figure  12:  Closed-Loop  Split  Model. 
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Figure  13:  Closed-Loop  Lumped  Model. 


dynamics  induced  by  the  control.  Consequently  designing  a  feedback  control  from  such  a  model  results  in 
very  ineffective  response  when  the  control  is  applied  to  the  system. 

Full  Order  Validation 

To  validate  the  effectiveness  of  the  reduced  order  control  obtained  via  the  split  basis  method,  we  utilize  a 
fixed-point  projection  algorithm16  to  incorporate  the  reduced  order  boundary  control  in  the  full  order  model. 
The  closed- loop  solution  of  the  full  order  system  is  shown  in  Figure  14.  As  seen  in  that  figure,  the  reduced 
order  control  effectively  drives  the  full  order  plant  to  the  target  profile.  The  full  order  discretized  model 
is  comprised  of  roughly  2,000  states.  The  reduced  model  obtained  via  the  split  basis  method  describes  21 
states.  As  a  result,  system  dimension  is  reduced  by  roughly  two  orders  of  magnitude  with  the  resulting 
reduced  order  control  being  quite  effective. 

Conclusions 


In  this  paper,  a  reduced  basis  construction  method  was  developed  allowing  for  separate  consideration 
of  baseline  and  actuated  dynamics  in  the  reduced  modeling  process.  A  prototype  initial  boundary  value 
problem,  governed  by  the  two-dimensional  Burgers  equation,  was  formulated  to  demonstrate  the  utility 
of  the  method.  A  weak  formulation  approach  was  used  to  develop  a  reduced  order  model  of  the  system 
with  boundary  control  appearing  explicitly  in  the  reduced  model.  When  actuated  energy  was  considered 
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Figure  14:  Closed-Loop  Response  with  Split  Method  Feedback  Control. 

separately,  much  better  agreement  was  seen  between  open-loop  solutions  of  the  reduced  and  full  order 
systems.  Separate  consideration  of  energy  induced  by  the  boundary  control  resulted  in  effective  feedback 
control  for  the  reduced  and  full  order  systems.  When  actuated  energy  was  not  explicitly  accounted  for  in 
the  reduced  modeling  process,  the  resulting  feedback  control  was  completely  ineffective  when  applied  to  the 
system. 

These  results  demonstrate  the  need  for  separate  consideration  of  baseline  and  actuated  energy  in  the 
reduced  modeling  process  when  the  resulting  model  is  to  be  used  for  feedback  control  law  design.  Basis 
construction  relying  on  an  energy  argument  applied  to  a  lumped  snapshot  ensemble  containing  baseline  and 
actuated  data  can  result  in  important  control  information  being  discarded  during  order  reduction.  This  is 
particularly  the  case  in  a  boundary  control  setting  where  ensemble  energy  is  typically  dominated  by  that  in 
the  baseline  data.  Feedback  controls  developed  from  the  resulting  reduced  model  are  likely  ineffective  when 
applied  to  the  system.  Separate  consideration  of  dynamics  induced  by  boundary  control  input  results  in 
reduced  order  controllers  that  are  much  more  effective  when  applied  to  the  reduced  and  full  order  systems. 
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